import math
from matplotlib import pyplot

gamma0=[0.1,0.3,0.5,1.0,2.0,5.0,10.0,15.0,20.0,30.0]
sigma0=[0.840309,2.2615,3.34583,4.77882,4.56534,3.21504,1.5779,0.931753,1.67977,16.3588]
delta=[1.824e-11,0.0016144,0.097718,0.263176,0.485487,0.528721,0.352905,0.138353,0.079131,0.372097]
omega=0.1

f = open("G1G2W30.0AMPL.dat","w")
for i in range(len(sigma0)):
    G1=100.0*sigma0[i]*math.cos(delta[i])/gamma0[i]
    G2=100.0*sigma0[i]*math.sin(delta[i])/gamma0[i]
    f.write(str(gamma0[i])+' '+str(G1)+' '+str(G2)+' '+str(sigma0[i])+' '+str(delta[i])+'\n')
    print "G': ",G1, "G'': ",G2
f.close()


#gamma0=[0.001,0.005,0.010,0.020,0.030,0.050,0.10,0.15,0.30,0.50,1.00,2.00]
#sigma0=[0.65711,0.332125,0.642324,1.15402,1.54863,2.00457,2.5053,2.86368,3.24615,2.81108,1.73702,4.29353]
#delta=[1e-11,1e-11,1e-11,0.104084,0.233288,0.440593,0.746536,0.867528,0.939347,1.00123,0.85841,0.710849]
#omega=0.1
#
#f = open("G1G2W10.0AMPL.dat","w")
#for i in range(len(sigma0)):
#    G1=sigma0[i]*math.cos(delta[i])/gamma0[i]
#    G2=sigma0[i]*math.sin(delta[i])/gamma0[i]
#    f.write(str(gamma0[i])+' '+str(G1)+' '+str(G2)+' '+str(sigma0[i])+' '+str(delta[i])+'\n')
#    print "G': ",G1, "G'': ",G2
#f.close()

gamma0=[0.1,0.5,1,2,5,10,15,30,40,50,80,100]
sigma0=[0.7008,2.50842,3.26897,4.03629,4.2139,3.14704,2.12561,10.7489,19.9096,25.1733,28.1566,19.994]
delta=[0.22845,0.473718,0.735112,0.941123,1.10058,1.11832,0.965732,0.944456,1.01967,1.0764,1.09221,1.15239]
omega=0.1

f = open("G1G2W10.0AMPL.dat","w")
for i in range(len(sigma0)):
    G1=100.0*sigma0[i]*math.cos(delta[i])/gamma0[i]
    G2=100.0*sigma0[i]*math.sin(delta[i])/gamma0[i]
    f.write(str(gamma0[i])+' '+str(G1)+' '+str(G2)+' '+str(sigma0[i])+' '+str(delta[i])+'\n')
    print "G': ",G1, "G'': ",G2
f.close()

gamma0=[0.01,0.02,0.05,0.1,0.2,0.5,0.7,1.0,2.0,5.0,7,10,15,20,30,50,100]
sigma0=[0.0526798,0.1097,0.264951,0.499303,0.875518,1.44475,1.62984,1.62568,1.85667,2.53464,2.92512,3.14699,3.52158,3.33532,3.65211,3.13998,3.54157]
delta=[0.388,0.4327,0.392697,0.41633,0.457377,0.611281,0.726921,0.925932,1.16584,1.36478,1.441,1.5354,1.57371,1.65698,1.61077,1.58688,1.44345]
omega=0.1

f = open("G1G2W2.0AMPL.dat","w")
for i in range(len(sigma0)):
    G1=100.0*sigma0[i]*math.cos(delta[i])/gamma0[i]
    G2=100.0*sigma0[i]*math.sin(delta[i])/gamma0[i]
    f.write(str(gamma0[i])+' '+str(G1)+' '+str(G2)+' '+str(sigma0[i])+' '+str(delta[i])+'\n')
    print "G': ",G1, "G'': ",G2
f.close()